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DISCOVERY 


SCIENTIFIC SOCIETY 


Investigating the lumped-heat- 
capacity assumption of thermal 


analysis in fibre mat composites 


Cornelius Ogbodo Anayo Agbo™ 


ABSTRACT 


The temperature gradient through the thickness of a composite material 
contributes to the thermal stresses experienced by the finally cured laminate. 
This paper seeks to evaluate the appropriate laminate thickness range where 
lumped-heat-capacity analysis can be applied. The conditions as evaluated 
made use of the Biot parameter. The material’s conductivity and the surface 
convection heat transfer coefficient coupled with the environmental 
conditions determine the borderline thickness for either using the lumped 
capacity analysis or the more detailed analysis. The variations in these 
parameters determine the temperature gradient through the thickness of the 
material. A numerical analysis carried out was validated using experimental 
tests. The results showed a negligible difference. The lumped-heat capacity 


assumption is therefore suitable in fibre mat composite analysis. 


Keywords: polymer composites, lumped-capacity, thermal properties, Biot 


parameter, thermal analysis 


1. INTRODUCTION 


The existence of spatial change in temperature in materials ensures that heat is 
conducted from one point to another within the material and the exchange of 
heat with the ambient. The temperature distribution in a laminate during cure 
influences the resultant residual stresses in the component and therefore 
require careful investigation and control. The temperature distribution in the 
laminate would depend on the material's thermal conductivity and the heat- 
transfer conditions from the laminate surface to the surrounding fluid. 
Reasonably uniform temperature distribution in the laminate would occur if 
the resistance to heat transfer by conduction were very small compared to the 
convection resistance at the surface so that the major temperature gradient 
would occur through the fluid layer at the surface. The lumped-heat-capacity 
analysis, then, assumes that the magnitude of the material’s external resistance 
is very much greater than the material’s internal resistance (Holman, 2001). 
Composite structures exhibit inherent processing defects some of which can 
be detrimental to the mechanical performance or constitute a nuisance than a 
significant problem; these defects include voids, delamination, residual stress- 
induced cracking, resin starvation, resin-rich pockets, damaged fibres, fibre- 
matrix debonding, thermal decomposition, and under cured and over cured 
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regions. The temperature field inside the fibre composite is very complex and directly affects the fusion quality between plies. There 
exists a characteristic nonlinear process in the internal temperature of the matrix as its temperature is raised during cure at it affects 
the quality of bonding of the composite layers and the component products (Cao et al., 2020). A major problem that could arise 
from the application of composite laminates is perhaps that of cracking due to large temperature differentials. 

Polymer composites generally are processed at increased temperatures, while thermoset matrix materials require elevated 
temperatures to cure, and thermoplastic matrix materials require elevated temperatures for the melting required for consolidation 
of layers. After a while at the elevated temperature state, the composite temperature is reduced, eventually reaching ambient 
temperature. Though the composite is soft and stress-free at the elevated temperature, at some point in the cooling process the 
material hardens to the point that the elastic and thermal expansion properties of the composite are identifiable. As the temperature 
of the composite goes below thermally induced stresses begin to develop (Hyer, 1998). Ultimately when the glassy state is reached, 
there remains a significant state of permanent stress, the residual stress state (Christensen, 1982). In the general case of heat flow, 
the temperature is not a straight-line function of the dimensions of a body. Since each element of a body is constrained by the other 
elements, it is not free to expand, and consequently, if the temperature gradient is not linear, internal stress will develop. Residual 
stresses are a system of stresses which can exist in a body when it is free from external forces and are generated by nonuniform 
plastic deformation (Dieter, 1988). The temperature gradient is likely to increase with the thickness of the laminate; therefore, 
greater thermal stress can be expected in thick laminates than in thin ones as the stress is directly proportional to the linear 
coefficient of thermal expansion, to the difference in temperature at the two faces of the laminate, and the modulus of elasticity 
(Agbo, 2016). 

After the cure cycle, deformations encountered during cure due to chemical and thermal changes coupled with the mould 
constraint are balanced internally in the component which results in residual stress that can cause low performance through 
warpage, shape distortions, matrix cracks and delamination within the manufactured part. An optimization process to predict 
accurately and measure accruing residual stresses and strain developments during the cure cycle constitute important issues 
needed to be addressed (Antonucci, 2006). 

Many polymeric composite laminates, like graphite/epoxy, undergo a substantial temperature change as they cool down after 
cure. The thermal excursions induce significant residual stress into the composite laminate because of the thermal expansion 
mismatch between the fibre and the matrix. The residual stress dependence on temperature history is attributable to two 
contradictory roles of temperature. While at elevated temperature stress relaxation occurs, upon cooling against geometric 
constraints, temperature acts as a stress-inducing agent (Weitsman, 1979). Linear elastic stress analysis predicts that the thermal 
stress, which develops within a composite material due to typical temperature excursions undergone by the structure, may exhaust 
the strength in various laminae even before external loads are applied (Lee and Weitsman, 1994). A major proportion of the thermal 
stresses is developed in the composite during cooling from cure temperature down to room temperature. This cool down stage can 
be considered as a part of the component manufacturing sequence and therefore lends itself to careful control. The cool-down path 
details appear to influence the magnitude of the residual stresses more noticeably than the duration of the cooling period. Thermal 
excursions induce significant residual thermal stress in composite laminate structures as an unavoidable problem of the large 
difference in the thermal expansion coefficient of the reinforcement fibres and the polymeric matrix. The presence of temperature 
gradient across the laminate thickness during cooling contributes to the development of residual stress in composites. This is due to 
the thermo-mechanical properties of the composites on temperature. As different sections of the composite laminate experience 
varying temperature histories, the deformations accruing therefrom will also differ. An inconsistency in the deformation response 
occurring in the laminate naturally leads to the development of residual stresses. Residual stresses may cause significant damages 
in composite structures, through void formation or fibre buckling during solidification, initiation of transverse cracks and 
delamination which reduces strength with a deleterious consequence on the component performance (Sonmez, 2002). 

In many studies done by earlier researchers, the temperature was assumed to be uniform across the plate thickness at any given 
time and so neglecting the effects of temperature gradient induced during cooling (Sonmez, 2002). Li et al. (2013) developed a 
model based on finite element analysis to predict residual stress occurring in thermoplastic composites. The values obtained 
indicated that both the sequence of the laminate layup and the distribution of crystallinity are affected differently by the cooling 
rate. Certain minimization of residual stress due to cure may result by understanding the cool-down path after cure. Although the 
residual stresses reduction in composites is not dramatic, they serve to indicate that additional benefits may accrue if a more 
scientific approach is taken toward the entire production cure process (Kibler, 1980). This paper, therefore, hoped to work towards 
predicting better temperature distribution criteria for a reduced residual thermal stress analysis. 
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2. THERMAL RESPONSE ANALYSIS 


Even though lumped heat capacity systems are considered uniform in temperature, these are idealized condition because a 
temperature gradient must exist within the material if heat can be conducted into or out of the body, in cooling a material within an 
environment, the spatial temperature distribution in the material would depend on the material’s thermal conductivity and the 
overall heat transfer conditions from the surface. A reasonably uniform temperature distribution exists when the resistance to 
conduction heat transfer was small compared to the convection heat transfer resistance from the surface so that the major 
temperature gradient would occur across the fluid layer at the surface. The lumped heat capacity analysis, therefore, presumes the 
internal conductive resistance of material negligible compared to external convective resistance. In that case, the convection heat 
loss from the body is evidenced by a decrease in the internal energy of the body (Holman, 2001). The classical lumped parameter 
analysis applies to various thermal engineering systems with low to moderate temperature gradient where the evaluated Biot 
number is below 0.1 (An and Su, 2015). A developed model of lumped capacitance was shown as a good replacement to the 
computationally complex CFD model in predicting the phenomenon associated with the transients in heat transfer (Tegenaw et al., 
2019). 


Heat Transfer Cooling Analysis 

Transient cooling analysis of the heat transfer is required to obtain the temperature distribution within the laminate and the 
temperature gradient through the laminate thickness during cooling. Necessary assumptions and boundary conditions were made 
to enable feasible diagnostics. The in-plane dimensions of the laminate were assumed to be very large compared to the thickness. 
Consequently, heat transfer through the edges has been neglected in comparison to the top or bottom surfaces of the laminate. The 
in-plane temperature distribution is accordingly assumed to be uniform. Thermal analysis of one-dimensional heat transfer is, 
therefore, considered and expected to yield acceptable results. 

The laminate’s upper surface is cooled by the convective heat transfer by the surrounding air, while the lower surface in contact 
with the mould is cooled by conduction heat transfer across the galvanized steel mould (see Fig. 1). The mould is in turn cooled by 
convection through the surrounding air. The conduction through the air film and radiation are all lumped into the convective 
overall heat transfer. The steel mould thickness is thin and highly conductive compared to composite laminate, hence the 
temperature at the laminate bottom is assumed to be the same as at the mould bottom. 


Laminate 
<< 
EE - mould 


Fig. 1 Laminate in mould undergoing cure process 


To establish an empirical value for the coefficient of heat transfer to be used, Nusselt number is sought as follows; 


The top surface of heated laminate as given by McAdams (1954): 


Nu, = 0.54(Ra)!/4 if 104 <Grr. Pri < 107 

= 0.15(Ra)!/3 if 107 < Grr. Pri< 10" (1) 
The bottom surface of heated laminate as given by McAdams (1954): 
Nur = 0.27(Ra)*/4 if 105 < Gr. Pre < 10% (2) 
Where Nu = Nusselt number 
Gr = gRL3AT/v? = Grashof number 


— ~ = Prandtl number 
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Ra = Raleigh number = Gr . Pr 


g = acceleration of gravity, v = kinematic viscosity, a = thermal diffusivity, 8 = 1/T is the volume coefficient of expansion, AT = 
temperature difference between the wall and the ambient, L= characteristic dimension, Q = density, t = time constant, 


The fluid properties are determined at the mean film temperature, 


= (To+T oo) 
2 


Ty where f indicates evaluation at film conditions 
Simplified equations for the heat transfer coefficients from surfaces to surrounding air at atmospheric pressure and moderate 


temperatures can also be approximated as follows (Holman, 2001): 


heated laminate facing upward, = h = 1.32(AT/L)1/4 laminar flow 104 < GrePre< 10° 
h = 1.52(AT)¥/3 turbulent flow  GréPri> 10° 
heated laminate facing downward h = 0.59(AT/L)1/4 both turbulent and laminar 


where h= heat transfer coefficient W/m2K 

AT = Ty —Two »°C 

L = vertical or horizontal dimension, m 
The laminate was kept horizontal on the mould to avoid resin overflow and/or work sliding down. The governing heat balance 
equation is stated thus: 

q = hA(T — T,.) = —cpVdT/dt 

where A is surface area, V = volume 
The initial condition 

T=Tp attime t=0 
So that the solution gives 


{7 > To) _ .-cmayecyyt 
(To _ Too) 


Where To is the temperature of the surrounding 

This analysis is expected to yield reasonable estimates within about 5% when the following condition is met. 
h(V/A)/k < 0.1 

where k= thermal conductivity. 

The dimensionless group is known as 

Biot number = hs/k = Bi 


S=V/A isa characteristic length dimension. 
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A representative element is shown in Fig. 2 to analyse the process of heat transfer within the composite. The Fouriers equation is 


established within the principle of the heat balance. 


qztdz 


Ok Qxt+dx 
Xx 
gen 
qytdy 
y 
Fig. 2 Heat transfer element from a laminate 
dT 
q = —kAS (3) 


Nevertheless, the general three- dimensional heat conduction equation is 


er eT aT gq 10T 


ie age oe ae o 


where the quantity a = = is called the thermal diffusivity of the laminate. 
For the one-dimensional case with constant properties, the differential equation for the temperature distribution 


T(z, t) in equation (4) is reduced to 


at = ade (5) 
the initial condition is T(z,0) = T; 

the convection boundary condition 

hA(Tup — T)z=9 = KASS at z=0 

The explicit forward difference equations are employed to predict the temperature distribution in the laminate as a discrete function 
of space and time. 

For the internal nodes, we have 
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: 2aAt 
TRY = 2 (Ther + Thea) + [1 as] me ©) 


ey 


and if the time and distance increments are chosen so that £ > 2, stability is ensured. 


The temperature of node m after the time increment is given as the arithmetic average of the two adjacent nodal temperatures at the 
beginning of the time increment. 

This restriction automatically limits our choice of At, once Az is established. 

On the convective surface of the composite laminate, the total energy conducted and convected into the node equals the increase in 


the internal energy at the node. Thus 


TP. —TP Az TPtt_TP 
kAy 2+ + hAy(T,, —T?) = pc — Ay 7 
sy em m) = pe Ay —— (7) 
aAt (_ hAz , [G2 hAz 
Tet = 2—-Ty +27? =f 217, 8 
me (Az)? k - at k (8) 


‘ : F ‘ Az)? . : 
To ensure convergence in the numerical solution, all selections of the parameter — must be restricted according to 


(Az)? naz 
mat = “GD 


In terms of resistance — capacity formulation (6) and (8) takes the form: 


At Tt) -T 
TP = at). J) (9) 


Where C; = p;c; AV; 
And q; = 0 for no heat generation 
Re = Ras = ~ inside the composite 


RSS at the surface of the laminate 


Table 1. Properties of air at atmospheric pressure 


TK Q Cp ut x10°, v x10° k a x104 Pr 
kg/m$ kJ/kg.K kg/m.s m2/s w/m.K m2/s 

100 3.6010 1.0266 0.6924 1.923 0.009246 0.02501 0.770 
150 2.3675 1.0099 1.0283 4.343 0.013735 0.05745 0.753 
200 1.7684 1.0061 1.3289 7.490 0.01809 0.10165 0.739 
250 1.4128 1.0053 1.5990 11.31 0.02227 0.15675 0.722 
300 1.1774 1.0057 1.8462 15.69 0.02624 0.22160 0.708 
323 1.0949 1.0072 1.9514 18.02 0.02798 0.2569 0.703 
350 0.9980 1.0090 2.075 20.76 0.03003 0.2983 0.697 
400 0.8826 1.0140 2.286 25.90 0.03365 0.3760 0.689 
450 0.7833 1.0207 2.484 31.71 0.03707 0.4222 0.683 
500 0.7048 1.0295 2.671 37.90 0.04038 0.5564 0.680 


Curled from (Holman, 2001) 
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The air properties change with pressure and temperature. At atmospheric pressure, the air properties vary with temperature as 


in Table 1 and visually presented in Fig. 3 for trend observation. The desired properties of air at 323 K is obtained through 


interpolation. 


——p 


—t=  x105, 
=< v x106 
=k 
== a x104 
== Pr 


Temperature (K) 


Figure 3. Variation of air properties with temperature, (K) 


The applications of equation (9) to laminates of usual thicknesses of between 3 mm and 10 mm as obtained within the auto 
industry will provide a good evaluation that may lead to general inferences. In the cases of laminate thicknesses of 3mm and 10 mm 
the temperature distribution across the thickness can be visualized from the following analysis shown in Tables 2, 3 and 4. The 
values in Table 2 is applied to determine the time increment appropriate for the calculations of the temperature distribution which 
is obtained as 3 secs. For the 3 mm thick laminate 4 nodes were chosen (see Table 3) while for 10 mm thick laminate, 11 nodes were 
chosen (see Table 4). 


Table 2. The determination of appropriate time increment for the 10 mm thick laminate 


Node Y/R) Ci Sam s 

1 208.912 627.48 3.00356131 
2 404 1254.96 3.10633663 
3 404 1254.96 3.10633663 
4 404 1254.96 3.10633663 
i) 404 1254.96 3.10633663 
6 404 1254.96 3.10633663 
7 404 1254.96 3.10633663 
8 404 1254.96 3.10633663 
9 404 1254.96 3.10633663 
10 404 1254.96 3.10633663 
11 205.456 627.48 3.05408457 
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Table 3. Numerical computations for the 3-ply 3 mm thickness laminate 


No. of time increment (At 
se), Ti T2 T3 Ta 

0 70 70 70 70 

1 68.6781 70.0000 70.0000 69.3391 
2 68.6766 69.3617 69.6808 69.3274 
3 68.0601 69.1850 69.3560 69.0189 
4 67.8887 68.7244 69.1107 68.6998 
5 67.4437 68.5074 68.7257 68.4571 
6 67.2336 68.0992 68.4906 68.0811 
7 66.8391 67.8702 68.1039 67.8474 
8 66.6175 67.4851 67.8672 67.4697 
9 66.2454 67.2506 67.4908 67.2344 
10 66.0185 66.8812 67.2510 66.8668 
100 47.2539 47.6987 47.8459 47.6920 
200 37.6239 37.8206 37.8855 37.8176 
500 30.6577 30.6747 30.6803 30.6745 
1000 30.0111 30.0114 30.0115 30.0114 
1200 30.0022 30.0022 30.0022 30.0022 


Table 4. Numerical computations for the 10-ply 10 mm thick laminate 


No. of time 

increment Ti T2 Ts T4 Ts To T7 Ts To Tio Tu 
(At =3 sec). 

0 70.00 | 70.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 
1 68.68 | 70.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 70.00 69.34 
2 68.68 | 69.36 70.00 70.00 70.00 70.00 70.00 70.00 70.00 69.68 69.33 
3 68.06 | 69.34 69.69 70.00 70.00 70.00 70.00 70.00 69.85 69.66 69.02 
4 68.04 | 68.89 69.67 69.85 70.00 70.00 70.00 69.93 69.83 69.44 69.00 
5 67.61 | 68.86 69.38 69.84 69.93 70.00 69.96 69.92 69.69 69.42 68.78 
6 67.57 | 68.51 69.35 69.66 69.92 69.95 69.96 69.83 69.67 69.24 68.75 
7 67.23 | 68.46 69.09 69.63 69.81 69.94 69.89 69.81 69.54 69.21 68.58 
8 67.19 | 68.17 69.05 69.46 69.79 69.85 69.88 69.72 69.51 69.07 68.55 
9 66.91 | 68.12 68.82 69.42 69.66 69.83 69.79 69.70 69.40 69.03 68.41 
10 66.86 | 67.88 68.77 69.25 69.63 69.73 69.76 69.60 69.37 68.91 68.38 
100 59.13 | 60.06 60.84 61.46 61.94 62.26 62.42 62.43 62.28 61.97 61.51 
200 53.03 | 53.76 54.37 54.87 55.25 55.50 55.63 55.64 55.52 55.28 54.91 
500 41.38 | 41.74 42.04 42.29 42.48 42.60 42.67 42.67 42.61 42.49 42.31 
1000 33.51 | 33.63 33.72 33.80 33.85 33.89 33.91 33.91 33.90 33.86 33.80 
2000 30.34 | 30.35 30.35 30.36 30.37 30.37 30.37 30.37 30.37 30.37 30.36 
3000 30.03 | 30.03 30.03 30.03 30.04 30.04 30.04 30.04 30.04 30.04 30.03 
3600 30.01 | 30.01 30.01 30.01 30.01 30.01 30.01 30.01 30.01 30.01 30.01 
4000 30.00 | 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 
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2. EXPERIMENT VALIDATION 

2.1. Materials 

E-glass chopped strand mat obtained as a continuous flat sheet serves to reinforce the polymer matrix. The fibres are evenly and 
randomly distributed and are held together in a powdered binder to obtain a mat area density of 450 g/m. It is quasi-isotropic in 
the in-plane directions. The polymer matrix is of the ortho-phthalic unsaturated polyester resin composition. The catalyst used is 
methyl] ethyl ketone peroxide, while the cobalt derivative served as the accelerator. Galvanized mild steel plate mould; oven; remote 
thermometer; scissors; brushes; knife; bowels; stirrer; metering syringes are part of the tools employed for the experiment. 


2.2. Methods 

The wet lay-up method was adopted during the laminate samples fabrication to ensure close monitoring and simplicity. Four pieces 
of mild steel sheets of 350mm x 220mm x 2mm were cut from a mild steel sheet and served as a mould. The laminate was prepared 
in two sizes, one has three plies of 320 X 150 mm and another has ten plies of 100 x 80 mm. The probe of a remote thermometer was 
securely attached to the laminate surface to monitor the curing temperature. Readings were obtained as the temperature rises 
during cure and falls after cure. The laminate was removed from mould 24 hrs later to avoid springing and bending that could arise 
if removed immediately after hardening. Three moulds of neat resin were also cured at different temperature regimes to appreciate 
its effects. 


3. RESULTS AND DISCUSSION 

From preliminary investigation, the maximum exothermic temperature reached by the laminate is 70°C at which cure takes place 
and also regarded as the free thermal stress state. The residual stresses are generated upon cooling to the ambient or storage 
temperature of 30°C 


704+30 


= 50°C = 323K 


Hence the film temperature Ty = 


The air properties at 323K were obtained by linear interpolation as shown in Table 1, therefore: 


_ (9.8)(0.003096)(70 — 30)(.06)? 


Ra = GrPr = aan toe (0.703) = 567526.992 ~ 5.68 x 10° 


Hence equations (1) and (2) applies 


Nu = 0.54(Ra)'/* = 0.54(567526.992)'/4 = 14.8214 


_ AL 


N 
CFE 


_ Nu(k) _ 14.8214(0.02798) 


top = —— ani = 6.912W /m? °C 
Nu = 0.27(Ra)!/* = 0.27(567526.992)!/4 = 7.4107 


7.4107(0.02798) ie 
a 50s = 3.456W/m2°C 


The condition for lumped capacity thermal analysis: 
: ‘ hs 

Biot number Bi= = < 0.1 

haverage = (6.912 + 3.456)/2 = 5.184W/m?eC 


for a laminate 3mm thick 
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__ 6.912 x 0.0015 


=0.05<04 
Bt 0.2020 


The result indicates that for a 3 mm thickness laminate lumped capacity can be comfortably applied as the Biot number obtained is 
below the 0.1 critical value. 


For a thickness of 10mm 


__ 6.912 x 0.005 


2 =0.17>01 
: 0.2020 


The result from 10 mm thick laminate is however above the specified critical value. Nonetheless, the difference is not large enough 


to distort the prediction arbitrarily, also, considering a more liberal condition that if 


lumped heat capacity analysis yields reasonable estimates. 


2x0.2020 


k, 
Thus, for all t < a ie( 
6h 66.912 


= 0.00974 = 9.74mm ~ 10mm) lumped-capacity analysis applies. The value determined from the 
second criterion is approximately 10 mm justifying the adoption of the lumped capacity assumption as well. 


Figures 4 and 5 show the change in temperature with time on the laminate surface with 3 mm and 10 mm thicknesses 
respectively as obtained through numerical analysis and by experimental measurements. It was observed that the numerical 
computation under-predicted the temperatures just after the commencement of the cooling process in the 3 mm thickness laminate 
while over-predicted the temperatures in the 10 mm thickness laminates. This may be due partly to the effect of edge losses 
predominant in thick materials unaccounted for in the thicker laminate analysis. Also, the abrupt changes in environmental factors 
idealized in the numerical method may have contributed to the discrepancy observed. The error is less in the 3 mm thickness 
laminate with 15% maximum as opposed to about 28% in the 10 mm thickness laminate. This suggests that lumped parameter 
analysis gives better estimates with thin laminates. 


x Experimental 
measurement 


@ Numerical computation 
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Figure 4. Cool-down temperature versus time for the 3-ply 3 mm thickness laminate 
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x Experimental 
measurement 
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Figure 5. Cool-down temperature versus time for a 10-ply 10 mm thickness laminate 


Figure 6 visualizes the temperature distribution across the thickness through the 10 mm thick laminate after 300 seconds of 
cooling down. The trend was determined from the values obtained through numerical computations. It was found that a 
preponderance exists in the temperature peak location towards the steel mould. This may be because heat loss from the top is more 
when compared to the heat loss through the bottom. 


Temperature (°C) 


4 6 8 


Thickness location (mm) 


Figure 6. Temperature profile across the 10-ply 10 mm thick laminate after 300 seconds. 


The associated biaxial stress in the laminate surface due to temperature gradient across the thickness is then obtained as follows 
(Hyer, 1998). 


=( 1+v JE AT 
~\GQ-v)a+v) « 
= EaAT 

1-v ” 
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PP = AT 
0, =03, =0 


EadT 10 
Toy Ee (10) 


Equation (10) gives the residual thermal stress within the composite as it cools from the high cure temperature down to ambient. 

where E is Young’s modulus, a = the linear coefficient of thermal expansion, AT is the temperature gradient and v is the 
Poisson's ratio. For E = 8 GPa, v= 0.4, a = 30x 10-°m/m/K. Considering the temperature distribution shown in Fig. 6 with a 
maximum change in temperature of 3.3 °C across the thickness from the core at node 8, and the surface at node 1. The maximum 
residual thermal stress due to temperature gradient is obtained as 1.32MPa. 

The major effect of uncontrolled temperature gradient in laminated composites is the residual stresses development which may 
lead to cracks and failure of the component even before been used. For an illustration, the cooled neat resin cured with cracks is 
depicted in Fig. 7. It is observed that as the temperature gradient increases the larger the cracks developed indicating higher 


residual stress accumulation. 


Figure. 7 A sample of cracking effect of cure temperature on neat resin mouldings. 


4. CONCLUSION 


It is established that the lumped parameter analysis subsists if the random fibre laminate thickness, the conductivity and the 
convective heat transfer are coupled within the Biot parameter and criterion. The spatial temperature distribution in the fibre mat 
reinforced laminate, therefore, can be considered under a thin plate with lumped parameter analysis if the thickness is below 10 
mm. For the fact that appreciable quantity of composite parts produced for the automotive industry fall within the thicknesses of 
consideration, 3 —- 10 mm (which is from 3 to 10 plies of 450g/m? glass fibre) it is sufficient therefore to apply the lumped-capacity 
analysis. An important remark though is the fact concerning the cooling and convective heat transfer input values and observations 
which were determined under a natural convective environment. This is so because most laminations using the wet hand layup 


method as prevalent in the automotive industries are cured under a natural environment. 
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